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q ■ Abstract 

We introduce a new concept , Stochastic Soliton Lattice, as a random process generated 
\Q [ by a finite-gap potential of the Shroedinger operator. We study the basic properties of this 
stochastic process and consider its KdV evolution 

> - 

1 Introduction 

o . 

The finite-gap solutions of completely integrable systems are almost periodic functions [N], 
[L], and, therefore, posess their natural stochastic structure determined by intrinsic shift 
group (compact group with the Haar measure) [PF]. Ergodic properties of finite-gap poten- 
tials were first used by Flaschka, Forest, and McLaughlin [FFM] in constructing the theory 
of multiphase averaging. 

The aim of the present work is the consistent introduction of the basic notions of the finite- 
gap stochastics, namely, we introduce a concept of a Stochastic Soliton Lattice (SSL) and 
study its basic properties. The stochastic structure of <?-gap soliton lattice (SL)(this term has 
been introduced by Dubrovin and Novikov in [DN] for finite-gap potentials) is the uniform 
distribution on g- dimensional phase torus (the real section of the correspondent Jacobian). 
We also generalize the important notion of the slowly modulated SL to the stochsatic case. 
The slow modulations of SSL are described by the hydrodynamic-type Whitham equations 
[W], [FFM], [LL], [DN] as well as in the case of usual SL. 

Following Johnson and Moser [JM] we define the rotation number for the SSL. It is well 
known that the rotation number coincides with the 'density of states' per unit length v(X) 
which is one of the most important quantities in physics of disordered systems [LGP] . For the 
SL the the rotation number is equal to the quasimomentum [Dl] that leads to a considerable 
effectivization of some general results of the theory of almost periodic functions. 

Passage to the stochastic description is especially natural when one considers the wave 
systems with a big number of degrees of freedom. 
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The case of the SSL with g » 1 and v « 1 (weakly interacting rare solitons Poisson- 
ically distributed on the line) corresponds to the 'soliton turbulence', studied by Zakharov 
in [Z] invoking the kinetic description. The similar object is studied by a different method 
in the paper of Gurevich, Zybin, and El [GZE] presented in these Proceedings. 

We note also that the introduced uniform phase distribution agrees with the Whitham 
approach to the description of nonlinear dispersive waves [W] , [DN] where the principal term 
of the WKB - decomposition (g-gap SL) is determined up to an arbitrary phase shift. 

Calculation of various averaged values on the realization space of SSL is an important 
applied aspect of the proposed theory. We present explicit formulas for the one-point proba- 
bility function and effective expressions for two first moments in the g-gap SSL. The challenge 
is to find the correlation function (the second mixed moment). 



2 Basic Preliminaries 

The soliton lattice (SL) is a function (g-gap potential of the Schroedinger equation) 

« p (x;r) = C7(r)-2^ B lne[z(x)|S(r)], (1) 

x e R, r = {n, . . . , r 2g+ i}, n < r 2 < . . . < r 2g < r 2g+1 ■ 

Here 

e[z|B(r)] = Yj exp{7ri[2(m, z) + (m, Bm)]} , (2) 

m 

m = {m 1 ,...,m g } e Z 9 , z e C m 
is the Jacobi theta-function of the hyperelliptic Riemann surface of genus g 

2g+l 

V 2 = Il( A -^) = ^ 5+ i(r,A). 
J'=l 

with cuts along the bands. 

The Riemann matrix B(r) and the constant C(r) are expressed in terms of the basis 
holomorphic differentials ipj : 



Bij 



Here 

9-1 

i>i = J2 a 3k I =d^ , (3) 

fc=o y/R 2g+ i(r,\) 
and dependence Ojfc(r) is given by the normalization 

j> i>j = S jk , (4) 
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where the a-cycles surround the bands clockwise, and the (3- cycles are canonically conjugated 
to a's. 

The imaginary phases z(x) are given by the formula 

z = — 2ia 3 _ 1 x + d , (5) 

where d is the initial imaginary phase vector. 
Using the substitution 

. = ^By (6) 
we rewrite ( [1]) in the form with real phases y [FFM] 

u g (x\r) =u g (y 1 (x),...,y g (x)\T), (7) 
yj(x) = kjx + fj , fj(mod2Tc) , (8) 

k = -27riS -1 a fl _i . (9) 

Here k = k(r) = {ki, . . . , k g } is the wave number vector, and f = {/i, . . . , f g )} is the 
initial (angle) phase vector, — n < fj < n. 
Note that in © 

u g (Vl, •••,% + 2tt, . . . , y g \r) = u g (y u . . . , y j} . . . , y g \r) , (10) 

that is u g (x\r) is (7-quasiperiodic on x (from here and on we consider only incommensu- 
rate kj,j — 1, . . . ,g. As for any almost periodic function we have for u g (x\r) the Fourier 
representation 

u g (x\r) = '£c j e i ^ x+h '\ (11) 
3 

where Cj , lj , hj are real, — 7r < hj < ir. (Really, lj G M, where M is a frequency modulo 
[JM] and hj are the integer linear combinations of fj(mod2ir)). 

The KdV evolution of (0) is isospectral and is described by the linear motion of the 
phases on the Jacobian: 

u g (x,t\r) = u g {y x (x,t), . . .,y g (x,t)\r) , (12) 
yj(x, t) = kjx + cjjt + fj , (13) 

where 

23+1 

oj = oj(r) = {wi,...,Wg} = -47riS~ 1 (a s _i J2 r i + 2a a-2) ( 14 ) 

i=i 

is the frequency vector. 
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We call u g (x\r) in (|7|), (|j) the modulated SL if r = r(ex) , e <C 1 (the possibility of 
= g(ex) is included). The initial value problem for the KdV equation 

u t - 6uu x + u xxx = , (15) 



u(x,0) = u g (x\r(ex)) (16) 

has the solution 

u(x,t) — u g (x,t\r(ex,et)) , (17) 

where u g (x, t\r) is given by (|r^),(|T3l) and the evolution of modulation parameters r(X, T) 
where (X,T) = (ex, eb) is given by the Whitham equations [W] , [FFM] , [LL] , [DN] 

drrj + Vjtfdxrj = , j = 1, . . . , 2g + 1 , (18) 



r J (X,0) = r,(X). (19) 

Here g = g(X,T) , g(X, 0) = g(X) and the Whitham equations (|18|),(|T9|) must be supple- 
mented with the additional information about evolution of g [LL], [VI], [DVZ], [Dl]. 



3 Stochastic Soliton Lattices: Definitions and Basic 
Properties 

It is well known that any almost periodic function generates the stochastic stationary process 
(see [JM], [PF]). 

Definition 1. The stochastic process generated by SL u g (x\r) we call Stochastic Soliton 
Lattice (SSL) and denote as v g (x\r). 

The general construction of u g (x\r) admits a very simple and clear description. The 
realization set Q of it consists of functions (0), (§) where f G T 9 ; T 9 is (/-dimensional torus 
(— 7r , n} 9 . The probability measure d\i is the uniform (Lebesque) measure on the torus (on 
Borel sets on it). It corresponds to the description of v g {x\v) following from(0),(§): 

v g {x\x) = u g {. . . kjx + 4>j... |r) (20) 

where 0i, . . . ,0 9 are independent random values uniformly distributed on (— n that is 
= {0i, . . . , <j) g } is uniformly distributed on T 9 . As kj are incommensurate then i> g (x\r) is 
an ergodic process [CSF]. 

As u g (x\r) is the stationary process then it has the Stone - Kolmogorov spectral decom- 
position [I]; due to the ergodicity this decomposition has the form (cf. (p|)) 

^|r) = E^ + ^, (21) 

j 

where 9j are the uniformly distributed on (— 7r , tt] noncorrelated random values [PR]. 



4 



The well known formula (Bochner - Khintchin) [I] gives us the covariance function K(h) 
of the stationary process v g (x\r): 

K{h) = (O g (x\r) ■ O g (x\r)) = £ |c,|V^ . (22) 

3 

Here = £(x) — (£) is the centered process. 

Theorem 1. Consider the KdV equation fll5|) as the equation describing the evolution 
in the phase space of stationary processes. Let the initial data has the form of the SSL: 

u(x, 0) = v g {x\v) . (23) 

Then the solution of the KdV is 

u(x, t) = u g (x\r) = u(x, 0) , (24) 

that is u(x,0) = u g (x\r) is a stationary point. 

Proof . The evolution of realizations ([7]) is described by flT2| ) ,(|Hf). For any moment 
£ one can introduce the new 'initial phase' f* = ujjt + fj which is also uniformly distributed 
on the torus T 9 . Therefore, the KdV evolution changes neither realization set Q nor the 
probability measure. 

Remark. In the connection with the Theorem 1 what can one say about the closure (in 
some topology) of the set {v g (x\r),r £ R 9 , g £ Z}? 



Definition 2. We call u g (x\r) from (20) a modulated SSL if r = r(ex). 
The modulated SSL (in opposite to the nonmodulated one) does evolve according to the 
KdV. 

Theorem 2. Consider the KdV intial value problem in the sense of the Theorem 1 with 
the same form of the initial data (|23|) but with r = r(ex). Then the solution of it is 



u 



'x, t) = i/ g (x\r(ex, et)) , (25) 



where r(X,T) is described by 

Remark . In a more general case the modulation function r(X) can be random itself (the 
simplest example is almost periodic r(X)). 

4 Some Mean Values (Moments) for SSL 

If Q(f) is an arbitrary smooth function then the mean value Q(£) of the stochastic process 
£is ' 



<Q(0> = / dfj.Q(f) , feQ 
h 



(We recall that £(x) = {Q = {f(x)} , B, //}, and B is some a-algebra of measurable Borel 
sets of Q). 
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For u g (x\r) we have [FFM] 

1 71 n 1 L 

(QMx\*)) = d<j> g Q{v g {x\v)) =lim-| Q(^(a:|r)) . (26) 



— 7T — 7T 



Direct calculation using (|T]),(0) gives surprisingly simple formulas for the mean value [V3] 
and the mean square of 

D g (x\r) = u g (x\r)-C(r). (27) 



Namely, 



and k(r),u;(r) are given by , ([HI) . It should be noted that expressions (p8[ ) are obtained 
for the particular choice of the canonical basis of cycles (see (|])). However, namely this 
normalizatiton is preferrable for many considerations of finite-gap solutions (see for example 
[FFM] , [V2]). 

We still do not know the expression for more important value: the second mixed moment 
(P g (x\r)u g (x + h\r)). 



5 One - Point Generation Function for the Moments 

The function 

/(A) = (e iXu ) (29) 

gives the moments (u) , (u 2 ) , . . . as the coefficients in expansion of /(A) in powers of A. 
Obviously, the one-point probability function is connected with /(A) by a Fourier transform: 
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Ug 



w(u) = — I f(X)e- iXu d\ (30) 

271 J 

— oo 

To calculate f g (\) for the SSL we make use of the trace formula for the SL u = %(y)© 
[Lev], [FFM] 

-r 2g+1 -2j2(H + Vj), (31) 

3=1 

,j = _r^i±^ (32) 

is the central point of the band and the functions [ij obey the Dubrovin's ODE's [D2] and 
each Hj lives inside the j -th gap. Then the ensemble integral (|26|) can be represented as 

7T 7T 



where 



-7T —77 
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f2 T2g 



The Jacobian ^ appears in [FFM] and has the form 



det * 



7r 9 d/x r2 ' +1 fc-i ' 1 J 

det[ / . rf/x] 
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Substitution of (^4|)into(|33|) allows to represent the multiple integral in the form of the 
product 

r2I + 1 e -2U(M+^) u j-i 

det[ / — dp] 



f g (X) = ' 2 ' r 2l / R ^ • (35) 



/jjtr 
dp] 

r 2l « 



Taking the Fourier transform of we get the probabilty density w q (-u)(|3"0"|). (As a 

oo 

matter of fact, / w g (u)du = 1 as f g (0) = 1). 

— oo 

The next step is to calculate the two-point generating function f g (Xi, X2] h) 

f 9 (Xi, A 2 ; h) = ( e -^M*) e -iMu(x+h^ ^ (36) 

which gives the covariance function (^) as the term before A1A2 in the decomposition of 
3611 



6 Density of States (Rotation Number) for SSL 

Consider the Schroedinger equation with almost periodic potential q(x) 

L<j> = (-d 2 xx + q(x))(f) = \<p, xeR. (37) 

The potential q(x) has a very important characteristics, the rotation number, which is 
defined for real A as [JM] 

a(X) = lim — arg(0/(x, A) — i(p(x, A)) . (38) 

We will use also the density of states p(X) which is connected with the rotation number 
by simple relation 

p(X) = -a(X) (39) 
7 



and the integral density of states [LGP] , [PF] , 



N(X) = / p(X)d\. (40) 



If q(x) is a finite-gap potential , q(x) = u g (x\r), then [DI] 

a( A)^, (44) 
where dp(X) is the quasimomentum [FFM], [DN] normalized by 

<fdp(\) = 0, j = l,...,g. (42) 

ft 

As the value of a (A) is the same for any realization from the set Q from Def.l then it 
is the spectral characteristics of the whole SSL v{x\r). Due to the Theorem 1 the rotation 
number does not change under the KdV evolution. 

We introduce also the full integral density of states in the SSL v as 

00 9 9 

u = N(oo) = - f dp{\) = — £ i dp(X) = . (43) 

IT J ATT J ATT 

— oo J~ L Otj J— x 

where kj are the components of the wave number vector @. Here we have used the known 
relationship [FFM],[DN] (recall that the a - cycles surround the bands in our normalization 
(see Sec. 2)) 

j> dp(X) = kj . (44) 

ay 

One can see that the full integral density of states has here the natural meaning of the mean 

number of waves per unit length. 

Remark. 

As the SL (0), (||)is the quasiperiodic function both in x and t we can deduce that an 
analogous result is valid for the full temporal density of states as well, namely 

^ = ^E«i. (45) 
where Uj are the frequences (|I4|), that implies that the temporal rotation number in the SSL 



is 



«-(A) = ^, (46) 



where dq(X) is the quasienergy [FFM],[DN] ( § dq{X) = ujj). 
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7 Distribution of Random Initial Phases in SSL. 
Poisson and Normal Limits 



We introduce the linear phases 

I JS * (47) 

where lj (j = 1, . . . , g) are independent random values uniformly distributed on (— £-] 
each. 

Introduce also the random value £ 3 - = X(o,i)(h) wmcn is the number of hitting of lj into 
the fixed interval (0, 1) C R (we suppose that ir/kj > 1). The variable £j takes two values: 
1 and with the probabilities Pj(l) = l/2n and Pj(0) = qj = 1 — Pj = 1 — kj/2n . The 
generating function <fj(z) of ^ is 

¥?i(z) = (l-Pi)+zPi (48) 

The sum £^ = X) is the number of hitting of all linear phases into (0, 1). As are 

i=i 

independent, the generating function for £^ has the form 

= n = na + (z - i) Pi ) = na + . (49) 
j= i j= i j= i zn 

Consider now the important particular case g>l There are two subcases to consider. 



7.1 Finite full integral density of states 

Suppose the following scaling 

k 3 = 0{g- 1 ) . (50) 
Then taking the logarithm of ( |l9f) we obtain 

hK P V>{z) = (z-l)±- y tk j + 0(g- 1 ). 

Therefore 

=exp{(z-l)z/}. (51) 

It is well known that the right-hand part of Q5TD is the generating function for the Poisson 
distribution with the mean value v. That means that is a Poissonic random value with 
the full integral density of states as a mean value. 
Remark 1. 

The scaling fl50l ) arises in particular in the semiclassical asymptotics for periodic KdV po- 
tentials, potentials [VI] [WK]. 
Remark 2. 

If v 1 then we arrive at the case of the separate oscillations asymptotically close to solitons 
Poissonically distributed on the line. The linear phases lj in this case give the asymptotic 
positions of the "solitons". This case corresponds to Zakharov's 'soliton turbulence' [Z]. 
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7.2 Large full integral density of states 



Suppose now that the following conditions fulfil 

1. Yf j= i kj = u^> 1, 

2. kj = o(l). 

Then the value ^ is distributed asymptotically as g — > oo by a normal distribution. This 
follows from the Lyapunov Central Limit Theorem in Berry - Essen formulation [F]. We 
present here the corresponding estimates for the distribution function F g (x) for the normed 
sums 

9 9 

E & - E p* 

k=l k=l 



The estimate is 
where 



(Ep*9*) V2 

k=l 

\F g (x) - <&(s)| = 0{g-W) , (52) 
$(x) = i- j e-y 2 /*dy. 
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